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Abstract 

In the Schrodinger formulation of non-Hermitian quantum theories a positive-definite metric op- 
erator n = must be introduced in order to ensure their probabilistic interpretation. This 
operator also gives an equivalent Hermitian theory, by means of a similarity transformation. If, 
however, quantum mechanics is formulated in terms of functional integrals, we show that the Q 
operator makes only a subliminal appearance and is not needed for the calculation of expectation 
values. Instead, the relation to the Hermitian theory is encoded via the external source j(t). These 
points are illustrated and amplified for two non-Hermitian quantum theories: the Swanson model, 
a non-Hermitian transform of the simple harmonic oscillator, and the wrong-sign quartic oscillator, 
which has been shown to be equivalent to a conventional asymmetric quartic oscillator. 

PACS numbers: 11.30.Er, 03.65.Db, ll.10.Kk 
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I. INTRODUCTION 



The recent interest in non-Hermitian quantum mechanics stems from the seminal paper 
by Bender and Boettcher[l|], who showed numerically that Hamiltonians of the form 

H =\p 2 - g{ix) N (1) 

possessed a real, positive spectrum for N > 2. This reality was attributed to an unbroken 
PT symmetry of H and was subsequently proved analytically by Dorey et al.Q, exploiting 
the connection between ordinary differential equations and integrable models. 

The more general framework of pseudo-Hermiticity, encompassing PT symmetry, was 
subsequently formulated by Mostafazadeh[3|, whereby 

H ] = 7]Ht]-\ (2) 

in which rj is a positive-definite, Hermitian operator. Essentially the same idea had appeared 
previously under the name of "quasi-Hermiticity" by Scholtz et al. [Ij , and indeed goes back 
to much earlier work by Pauli[5j]. 

Reality of the eigenvalues is not sufficient for a viable non-Hermitian quantum theory: 
for a probabilistic physical interpretation we must also have a positive-definite Hilbert-space 
metric. Such a metric may be constructed, but it is not given a priori, and depends on 
the particular Hamiltonian H. In the context of PT symmetry a new operator, the C 
operator, was introduced in Ref. [6[, with the help of which a positive-definite, CPT inner 
product could be constructed. In a subsequent paper 0] a systematic way of constructing C 
in perturbation theory was developed, which was greatly facilitated by the introduction of 
the Q operator, defined by PC = e~® . The connection with the operator rj of Eq. ([2]) is 

77 = e~ Q . (3) 

This appears in the overlap (( , )), defined by 

«^» = (^e-^>, (4) 

and in the expression for the equivalent Hermitian Hamiltonian j3[: 

h = e~^ Q He^ Q , (5) 

showing that the two operators are related by a similarity transformation (rather than a 
unitary transformation). 

In most cases one must resort to perturbation theory or some other approximation in 
order to construct Q and h. However, in the two theories we study below they are known 
exactly. These are the Swanson[8| Hamiltonian, 



H = uja)a + aa 2 + (3a) 2 , (6) 

where a and a) are simple harmonic oscillator annihilation and creation operators for unit 
frequency and u, a and (3 are real, and the "wrong-sign" quartic oscillator, 

H = \pl + \m 2 z 2 - gz\ (7) 
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where we have written z, p z in place of x, p in recognition of the fact that for the energy 
eigenvalue equation to be well defined, it must be formulated on a contour in the lower- half 
complex z plane. 

The Q operator (or operators, because it is not unique) can be found exactly (lol. Q | for 
Eq. and the equivalent Hermitian Hamiltonian is simply a harmonic oscillator. In the 
case of Eq. fl7|) it can again be found exactly for a particular choice of contour |9j, leading to 
an equivalent Hermitian Hamiltonian which is a standard quartic oscillator, together with 



a linear term, as discovered earlier by Andrianov[12| and Buslaev and Grecchi[13 



Apart from the standard formulation of quantum mechanics in terms of operators x and 
p acting on Schrodinger wave functions we may also approach it from the point of view 
of functional integrals, which indeed is the natural language if a generalization to higher 
dimensions is to be attempted. A brief discussion of the Swanson model in this context was 
given in Ref. [3] and a rederivation of the results of Ref. || using functional methods was 
given in subsequently corrected in 16]. These papers were primarily concerned with the 
derivation of the equivalent Hermitian Hamiltonian h and did not address the question of 
overlaps or expectation values. Indeed, the quantity considered was the partition function, 
i.e. the vacuum-generating functional for j = 0, rather than the true functional Z[j], from 
which Green functions can be obtained by functional derivation. 

In the following two sections we will explore the full vacuum generating functional, finding 
that in this formalism the operator Q appears only fleetingly and is not involved in the 
calculation of overlaps or expectation values. It is rather the dependence of Z[j] on j that 
encodes the similarity transformation induced by Q. In passing, we also note the role of Q 
at the classical level as being a generator of canonical transformations. In the final section 
we discuss the implications of our results. 



II. Z[j] FOR THE SWANSON MODEL 



In terms of x and p the Hamiltonian of Eq. reads 

H = ax 2 + bp 2 + c{x, p}, (8) 

where a = ~(u + a + (3), b = \{oj — cx — (3) and c = \i(ot — /3). It is non-Hermitian for a ^ /3. 
The Euclidean phase-space functional integral for Z[j] corresponding to Eq. (|S]) is 

Z\j) = J [D<p][Dic] exp { - J dt[i(pir + atp 2 + bn 2 + 2apn - jcp]}, (9) 

where, to correspond to the notation in field theory, we have written (p for x and it for p. 

The different solutions for Q and h are quite simply obtained by completing the square 
in the integrand in different ways 1 



Below we have illustrated the two extreme cases where Q = Q{(fi) and Q = Q(ir). The more general case 
Q = Qi^i 77 ), in which both tp and 7r are shifted, does not add anything essential to the case Q — Q(ir). 
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A. Q = Q(ip) 

This corresponds to completing the square by shifting it, according to 



ap 2 + bir 2 + 2cpii = b(ir + -if) 2 + (a — ^")^ 2 = bit 2 + cup 2 , 



(10) 



giving 



[Dp][D-K] exp | — J dt[ip(TC — -p) + aif 2 + bir 2 — jp]} 
[Dip] [Dn] exp | — J dt[iip7r + dp 2 + bn 2 — j 



The Q(x) operator has previously been determined in the operator formalism 11 



as 



Q{x) 



a — (3 
uj — a — P 



x 



-i—x 
2b 



(11) 



(12) 



Its time derivative, written in terms of p, appears in the first line of Eq. ffTTl) . but then 
disappears in the following line, precisely because it is a total derivative. The momentum 
variable can now be integrated out 2 , to leave 



. 2 

Z[j] = I [Dp] exp { - J dt + dp 2 - jp] } 



(13) 



the standard form for the vacuum generating functional of the simple harmonic oscillator 
with frequency Q 2 = Adb. Note that by completing the square in this particular way, we have 
effectively changed 7r to ff = 7r + (c/b)p but left p unchanged. In the operator formalism 
this would be the result of the transformations 



P 

x 



P = e ^ Q{x) p e ~^ Q{x) =p+ -iQ'[ 



x) 



X 



e2 Q{x) x e \Q(x) _ x ^ 



(14) 



In terms of classical mechanics this is a canonical transformation, generated by 

F 2 (x,P) =xP-Uq(x) 

according to the standard equations 

v dF 2 dF 2 



OP 



dx 



(15) 



(16) 



It transforms the original non-Hermitian Hamiltonian H of Eq. (jSJ) into the equivalent Her- 
mitian Hamiltonian 



h = bP 2 + dX< 



(17) 



Here, and in subsequent integrations, we ignore overall multiplicative factors. 
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while the corresponding Lagrangians, i and L, differ precisely by the time derivative of Q. 

In this way of constructing h the field ip remains an observable in the original theory, and 
Green functions of <p are obtained in the standard way by functional differentiation with 
respect to j(t). Since the single- variable functional integral for Z[j] obtained from Eq. is 
identical that of Eq. (1131) . the Green functions of the Swanson Hamiltonian can be obtained 
in the functional formalism without reference to Q, in contrast with the corresponding 
calculations within the standard Schrodinger framework, where matrix elements have to be 
calculated with the metric 77 = e~ Q , as in Eq. (0J. 



B. Q = Q(tt) 

In this case, there will be a shift in (p between the original and the Hermitian Hamiltonian. 
For that reason, and to conform with the following section, we change the notation so that 
the field appearing in Eq. (Q becomes ip: 

Z[0] = J [Dtp][Dn] exp { - J dt[itpn + aip 2 + bn 2 + 2apn]}, (18) 

reserving ip for the shifted field appearing in the Hermitian version. Note that for the moment 
we have set j = in Eq. (fl8"j) . This is because we will ultimately wish to calculate Green 
functions of the observable ip, and hence will insert j in the latter form and then transform 
back to determine the j dependence in Eq. ([TBI . Completing the square according to 



aip 2 + biT 2 + 2cifjir = a(ip + -vr) 2 + (b )tt 2 

a a 



= aip 2 + &7T 2 , (19) 

we obtain 

Z[0] = J[D(p][Lhr]exp{- J dt[i (<p - ^-7r)n + a<p 2 + 6vr 2 ]} 

= / [Dip] [Dtt] exp {- J dt[iipn + aip 2 + 6vr 2 ]}, (20) 

after discarding the derivative of Q(tc) = —m 2 {c/2a) appearing in the first line of Eq. ( 1201 ). 
Performing the integration over it we obtain 

Z[0] = J [Dip] exp { - J dt + aip 2 ] }, (21) 

which is a rescaled version of Eq. f|T3|) with j = 0. The corresponding operator transforma- 
tions would be 

p -> P = e^ Q{p) p e -^ Q(p) = p 

x — > X = e* Q(p) x e'^ Qip) = x- -iQ'(p), (22) 
while the corresponding transformation in classical mechanics would be generated by 

F 2 (x,P)=xP-±iQ(P). (23) 
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In order to obtain Green functions of ip we now restore j to Eq. (l2Tj) : 



dt 



-~+a(p -jip 
4b 



}■ 



Zbl = / [Df] exp { - y , 
or equivalently in Eq. (1201) : 

= /" [Dtp] [Die] exp { - J dt[i(pir + atp 2 + bn 2 - jip]}, 



(24) 



(25) 



and transform back to the original field variable ip by substituting (p — ip + (c/a)ir. The 
resulting expression in the exponent: 

| ] = iipn + bn 2 + aip 2 + 2cipir — jip — j (c/a)ii 

can be written in the form 



7T + ^ ( II ' + 2cl ' - -j 



c 



1 ', o ~ t 2 ^ ■ / ■ / 

H ib + ayj iw H 



which gives 

Z[j] = j [Dip] exp{-| 



dt 



1; 2 1„ 2 , 2 O 2 J , JC j ; C 2 J 2 

+ -Q 2 V> h -=ip — 

2^ 2 ^ 2a a^/W 4a 2 b 



(26) 



(27) 



after performing the 7r integration and rescaling ip by ?/> — > ipy2b. 

Let us now verify that we obtain the correct expressions for some low-order Green func- 
tions in (p by functional differentiation of Eq. (12"7|) . A single differentiation 8/8j(t) gives 



1 5Z 

zJJ 



j=0 



2&\2a o7' 



(28) 



which is clearly zero, the correct result for ((p). 

For a somewhat less trivial check we perform a double differentiation 8 2 /5ji5j 2 , using an 
obvious subscript notation for j at the two times ti, t 2 . Thus 



1 5 2 Z 



Z 5j 1 5j 2 



j=0 



'n 2 , ic • 

—1p 2 1p 2 

2a a 



+ 



2a?b 



S(t!-t 2 ), 



(29) 



This is simply evaluated by noting that for the harmonic oscillator Lagrangian appearing in 
Eq. (1271) the expectation values are 



The result is 



1 5 2 Z 
Z 5j 1 5j 2 



1 



,-flltj-tal 



2Q 

L(t 2 _ t^e-ni*!-*.! 



(30) 



i=o 



~2 e 



— < 

4a 



-fi|ti-t 2 



+ - t 2 ) 



-fi|tl-*2| 



(31) 
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after using the relation c 2 + Q 2 /A = ab. On the other hand, from Eq. ff24l) we expect 

(vm) = ^~ nltl ~ t2l > ( 32 ) 

which is indeed the same because Q 2 = Aab. 

The lesson we draw from these simple calculations is that the Green functions of the 
observable field (p appearing in the Lagrangian form of the equivalent Hermitian Hamiltonian 
h can be calculated using the functional integral arising from the non-Hermitian Hamiltonian 
H without any reference to the Q operator, which is needed to calculate expectation values 
in the Schrodinger formulation of H. Instead, the information about the transformation 
induced by Q is encoded in the functional dependence of Z[j] on j. 



III. Z[j] FOR V(z) = -gz 4 

As already mentioned, the Hamiltonian of Eq. (J7j) gives rise to a well-defined eigenvalue 
problem only if it is defined on a suitable contour in the lower-half complex z plane. The 
essential requirement on this contour is that it must lie asymptotically within the Stokes 
wedges [ll], which in this case extend from the real axis down to an angle of tt/3. Along 
the centre of the wedges the wave-function behaves purely exponentially, changing to purely 
oscillatory at the edges. The particular contour chosen in Ref. (9[, 

z = -2iy/(l + ix), (33) 

happens to go to infinity at an angle of — 7r/4, i.e. not down the centre of the wedge, but 
nonetheless it has some very special properties that enable an exact evaluation of the Q 
operator and the equivalent Hermitian Hamiltonian. The results are that 

Q = -^ + 2 P , (34) 



where a = lQg, and 



lba 4 



Notice that this has the unusual feature that p appears to the fourth power. However, by a 
Fourier transform h reduces to a standard, asymmetric quartic anharmonic oscillator. Note 
also that because Q = Q(p) it induces the transformations of Eq. (122]) . whereby 

did , . 

P = p = - l T x = -zTz- (36) 

The result for h has been reproduced in the functional integral formalism 15, [HI], but again 
the quantity considered was Z[0], not the full vacuum generating functional Z[j]. In the 
present section we review the functional integral derivation, and include an external source 
in the Hermitian formulation, tracing it back to the non-Hermitian functional integral. 
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A. Z[0] 

We start with the single functional integral 

Z = J[Dip] exp j~y" dt 

defined, as discussed above, on the curve C given by 

ip = -2iy/(l + i<p) 



^ip 2 + ^m 2 ip 2 - gip 4 



(37) 



(38) 



where <p is real. Making the change of variable from ip to ip, and including the additional 
effective potential induced by this non-linear change (9I. we obtain 



Z[0] 



[Dip] 



Dety/(1 + icp) 



exp 



dt 



1 2 



2 1 + i<p 32 1 + iip 



- 2m 2 (l + i(p) - a(l + i<pf 



Now we represent the functional determinant by a functional integral over 7r: 

IN 2' 



1 



DeV(l+i^>) 
thus producing the phase-space functional integral 

Z[0] = J[D(p)[Dir)exp{- J dt 



[Dn]e W {- ldt-(l + iip) l*-T^i£ 



(39) 



(40) 



-(1 + i(p)ir 2 — iipTr 2m 2 (l + iip) — a(l + i(p) 2 

2 4 



(41) 



Here the expression in square brackets in the exponent can be written as 



a 



(¥>-*') + — (-( 
Ky J 2aV2 v 



7r 2 - 4m 4 ) + 7T 



7T 



2\2 



7r (7r 2 - 4m 2 ) 



4a 4 16a 

So after rescaling ip — > Lp/\ / 2a, tc — > 7r\/2a and performing the integration we obtain 

^2 



(42) 



Z[0) = J[Dn]e W l-Jdt y - /| 



2m 2 \ 2 



/a a / 

' — 7T H 7T 

4 I a 



(43) 



in which (with p <-> 7rv2a) we can recognize the Lagrangian corresponding to the Hermitian 
Hamiltonian h of Eq. (1351) . Note again that in the transition to this last equation Q put in 
a brief appearance in the form of its derivative Q(7r), which was immediately dropped. 



B. Z\j) 

Now let us couple tt to an external source j by the addition of a term —jit to [ ] and 
work backwards to see how j appears in the the original functional integral. In Eqs. fj42|) 
and ( l4TT) the term 7r/4 is replaced by 



7T 7T 

4^4' 



1 + 4]), 



(44) 
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where j = j / \/2a. Then the 7r-dependent terms in the new version of Eq. (|4ip can be 
written as 



-(l + ^)vr 2 



UfTC 



4(l + 4i) 



1 / ^+i(l + 4j) 

— (1 + Hp) 7T 

2 V y! \ l + i(p 



^ 2 



1 



/ 2 1 + i(p 32 1 + zip 
So in the [ ] of Eq. ( 1391) we will have the additional terms 

Aiifj +J + 2] 2 



jupj + 3 + iy 
4(1 + i<p) 



(45) 



4(1 + i<p) 

Writing these in terms of ip using 4(1 + icp) = —tp 2 , so that itp 



(46) 



— ^tpip, we finally get 



Z[j] = / [DiP}exp 



c 



dt 



l(*P 2 + m^) - -] (2^ - ±) + 2 4 



A functional differentiation with respect to j gives 



1 5Z 



(vr) 



1 



3=0 



(47) 



(48) 



Unfortunately, although we are assured that H and h are equivalent, with identical energy 
spectra, we are unable to solve either exactly. Likewise we are not able to evaluate (71) 
exactly in either theory. However, in the next subsection we will check Eq. (1481) up to first 
order in perturbation theory. 



C. Perturbation theory for (it) 
1. Hermitian formulation 

First we calculate (71*) using normal Schrodinger perturbation theory for the Hermitian 
Hamiltonian h of Eq. fl35|) . As a change of notation, let p = yV2a, so that x = p y /\^2a, 
and h now reads 

1 9 1 ( 9 2m 2 \ 2 fa , jnS 

h= 2 p ' + r{» -—) -m (49) 

In order to perform perturbation theory in the standard way we need to shift y according 
to y — > y + (3 so as to eliminate the linear term. This yields the cubic equation 

a/3 3 - 2m 2 (3 - J- = 0. (50) 

Although in principle it is possible to solve this equation exactly, it is extremely cumbersome 
to do so: instead we will make successive approximations in powers of a. In fact, the lowest- 
order approximation is that (3 is large, of order 1/y/a. We thus write (3 = X/y/a, to obtain 
the equation 

X 3 -2m 2 X = ^=, (51) 
2V2 
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whose solutions are approximately A = ±my/2, — a/(4m 2 \/2). The overall minimum is at 
the first of these, A = m\/2, and the lowest-order mass term is ^M 2 y 2 , where M = 2m. 
Thus the zeroth-order result for (n) is given by the classical position of the minimum, i.e. 

<7T>=/? = -^L. (52) 



In order to calculate the first-order result we need to do two things: first calculate the posi- 
tion of the classical minimum to higher precision, and then calculate quantum corrections. 
Calculating the minimum to next order we find that 

where we recall that a = 16g. The shifted Hamiltonian now reads 

h=\p 2 y + \M 2 y 2 + mv^ (l + ^)v 3 + \ay\ (54) 

which gives a quantum contribution to (y) from a tadpole diagram arising from the cubic 
coupling. Its contribution is 

-M\t-t'\ l Sm^to 



A(y) = -3mV2a / dt' — — = i— 55 

xy/ J 2M 2M 2M 3 v ; 

Combining this with the classical shift of Eq. fl53|) we finally obtain 



\f2a V m J 



2. Functional integral formulation 

Now we wish to rederive these results from the functional integral of Eq. fT4T|) and the 
expression for (tt) of Eq. (|48l) . In the latter equation a considerable simplification can 
be made by noting that ip/ip is a perfect time derivative, and hence its expectation value 
vanishes, since there is no preferred origin of time. The remaining term, l/ip 2 , might seem 
singular, but we recall that the contour C stays away from the origin. In fact C is not 
suitable for performing perturbation theory, because it is right on the edge of the Stokes 
wedge for the harmonic oscillator. But, and this is the great advantage of the functional 
integral approach compared with the Schrodinger approach to the non-Hermitian problem, 
the contour may be analytically continued 3 to a contour compatible with the asymptotic 
requirements of both the harmonic oscillator and negative quartic oscillator potentials and 
still avoiding the origin. 



3 This is not possible in the Schrodinger approach, because the calculation of expectation values also involves 
the complex conjugate of the wave-function. 
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The calculation of (it) to lowest order can be achieved by the use of functional integration 
by parts. We start from the identity 







8 



-i J dtdt'il>(t)G- 1 (t-t')iP(t') 



5(tt - 1 2 ) 1 



dt'Go^ti-t'Wt') 



e 2 



(57) 



- 1 



| f dtdt'tP(t)Gg 1 (t-t')lP(t') 



where G (t — t') is the operator (—Of + m 2 )5(t — t'). Integrating this with respect to 
/ dtiGo(t2 — ti), we obtain 

= J [#] 

Since G (t — t') = e~ m ' i_ *''/(2m), this proves that 

/ 1 ' 



go(0) 



(59) 



^ 2 



-2m 



(60) 



and hence (7r) = 2m/v / 2a, in accordance with Eq. (1521) . This last equation, Eq. (150]) . might 
seem rather odd, but it must always be remembered that if) is complex. In order to obtain 
the first-order correction, we must keep the interaction term in the exponent. The analogue 
of Eq. (1571) is now 







-e 2 



| / dtdt' tpi^G' 1 {t-t')ip(t')+g J dtip 4 



S{t!-t 2 ) 1 



dt'G^ih-t'^^+Ag^y 



(61) 
(62) 



Now integrate / dtiG (t 2 — ti) to obtain 
= / [(ty] 



G ° (0) 1 + 4^ / dtiGofo - WW 



i.e. 



' 2m + -^r / d*ie- m|ta - tl| ^(*i) 3 



V^ 2 ) 2 



(63) 



(64) 



Given the fact that in the Hermitian sector we found that 2m — > 2m — 4g/m 2 , we need to 
show that 



(65) 



This is an unfamiliar kind of expectation value, but we can represent it in terms of more 
familiar quantities by using the integral representation 



4>(h) - ^ Jo 



(66) 
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In fact, precisely because our contour does not go through the origin, we may drop the is 
in the calculation of J, which now reads 



r , , r°° 00 (iX\ 2n+1 1 \ 

^-ij^e^l S ^L_(^>. (67 ) 

Thinking graphically in terms of contractions, these can occur in two ways: 

(a) («i)(^ 2 )(^ 2 r 

This occurs with coefficient 3(2n+ 1)!! and value e m ' 41 * 2 l/(2m) n+2 

(b) (^ 2 ) A (^M n - 1 

This occurs with coefficient 2n(2n + 1)!! and value e ~ 3,n l* 1_ * 2 l/(2m) n+2 

Thus 

j a = J* /°° d tl e- 2m ^ r dX \e-" 2 /^ = -L- , (68) 
4m z Jo Jo 2m z 

while 

h = l — r ^ie- 4m|tl " t21 r dX X 3 e~ x2 ^ = l — . (69) 

8m 6 Jo Jo 2m z 

Thus I = I a + If, — 1/ m 2 , in agreement with Eq. fl65l) . 

An alternative way of proceeding, which of course gives the same result, is to treat / 
in terms of a vacuum expectation value of a time-ordered product of operators and insert 
a complete set of states. The states involved are only the first and third excited states, 
and the required matrix elements of I/1P2 can be calculated from those of ip2 by writing 
((V^) 27 *) = {{'4 ) 2) 2r+l x (1/^2)) f° r r = 0, 1 and inserting the requisite number of intermediate 
states. 



IV. DISCUSSION 

The main thrust of this paper has been to illustrate how the metric, in the guise of 
the Q operator, makes only a fleeting appearance when the theory is formulated in terms 
of functional integrals rather than the operator formalism. Green functions are calculated 
as functional integrals in the normal way. The information about observables, namely the 
transforms of Hermitian operators in the equivalent Hermitian formulation, is carried in 
the non- Hermitian version by the functional dependence of Z[j] on the external current 
coupled to those observables in the latter formulation. In the case of the wrong-sign quartic 
oscillator, there is in fact a technical advantage in the functional formulation compared with 
the normal Schrodinger formulation, namely that the original contour can be distorted so 
as to make perturbative calculations of Green function possible. 

However, it must be admitted that the number of non-Hermitian theories, regarded as 
fundamental theories, that are interesting and tractable, is at the moment rather limited. 
One other theory that has been the subject of extensive investigation is the igx 3 potential. 
In this case one can stay on the real line, but the Q operator, and hence the equivalent 
Hermitian Hamiltonian h, is only known in low-order perturbation theory (IH. jToj ] . It has 



been shown that perturbative calculations of the energy levels are considerably easier in the 
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context of the non-Hermitian theory, regarded as a field theory in 1 dimension, than in the 
equivalent Hermitian theory, where the potential in h contains a number of momentum- 
dependent terms 20]. However, calculations of Green functions are likely to be difficult in 
both versions (2l|. and indeed by the methods of the present paper, because the transforma- 
tions induced by Q are complicated, with Q being a function of both x and p. 

An alternative way of viewing non-Hermitian theories is as effective theories, which hap- 
pen to be useful in deriving tractable results, the prime example being the non-unitary Dyson 
mapping used by Snyman and Geyer[22|] to convert the Hermitian Richardson Hamiltonian 
into an equivalent non-Hermitian boson-fermion Hamiltonian. 

Nevertheless, the — gz A potential is still of great interest as a fundamental theory because 
of its possible generalization to a 4-dimensional — g(p A theory. The reason is that the indi- 
cations in perturbation theory are that such a theory could be asymptotically free, as was 
already stressed in the original paper of Bender and BoettcherpJ, and spelt out in somewhat 
more detail in Ref. [13]. The first attempts at a generalization to higher dimensions were 



made in Ref. [151 . but it remains a challenging problem. 
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